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Abstract 

By the methods of multitype branching processes in random environ- 
ment counted by random characteristics we study the tail distribution of 
busy periods and some other characteristics of the branching type polling 
systems in which the service disciplines, input parameters and service time 
distributions are changing in a random manner. 
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1 Polling systems with service policies of branch- 
ing type 

We consider a polling system consisting of a single server and m stations with 
infinite-buffer queues indexed by z G {!,..., m} . Initially there are no customers 
in the system. When customers arrive to the system the server starts immedi- 
ately the service by visiting the stations in cyclic order (1— >2— >----^m— > 
1 — » ■ • • ) starting at station 1 according to a selected service policy (to be de- 
scribed later on) and with zero switchover times between queues. Later on the 
initial stage of services (1 ^ 2 ^ • • ■ ^ m) will be called the zero cycle. The 
subsequent routes of the server will be called the first cycle, the second cycle 
and so on. When the system is empty the server waits their arrival at a parking 
place R. Customers arrive to the queues in accordance with a point process 
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whose parameters are changing in a random manner each time when the server 
switches from station to station . 

To give a rigorous description of the arrival and service processes for the 
system in question we need some notions. Let s : = (si, . . . , Sm) G [0, 1]"* be a 
m-dimensional variable, 



Gj := (0, . . . , 1, . . . , 0) bo a m— dimensional vector with zero components except 
the i-th equal to 1, and 1 are m— dimensional vectors whose all components 
are equal to and 1, respectively. Let 



be the mixed probability generating function (m.p.g.f.'s) of (m + 1)— dimen- 
sional vectors {On, . . . , 9im', <Pi), where % are nonnegative integer- valued random 
variables and (f)i is a nonnegative random variable. Denote 



the respective vector-valued m.p.g.f. and let = {(f){s;X)} be the set of all 
such m.p.g.f. Assume that a probability measure P is specified on the natural 
cr- algebra of !F and let 



be a sequence selected in iid manner from in accordance with P, where 



In the present paper we investigate such polling systems whose arrival and 
service procedures of customers meet the following property. 

Branching property. When the server arrives to station i for the n-th time 
and find there, say, ki customers labelled 1, 2, fcj , then, during the course of 
the server's visit, the arrival of customers is arranged in such a way that after 
the end of each stage of service of customer j (the number of such stages may be 
more than one if the service discipline admits feedback) the queues at the system 
will be increased by a random population of customers {0ii{n,j), . . . , Oim{n,j)) , 
where 9ii{n,j) is the number of customers added to the Z-th station, and, in 
addition, a final product of size (j)i{n,j) > will be added to the system. It is 
assumed that the vectors {0ii{n,j), . . . , 6im{n,j)] (j)i{n,j)), j = 1, 2, ... , ki are 
iid and such that 




s^-,fcieNo = {0,1,2,...}, 




<^(s;A) := (</)W(s;A),...,<^W(s;A)) 



(/)o(s;A),^i(s;A),(/)2(s;A),... 
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Note that (j)ii'{s-,X) is a random m.p.g.f. and, therefore, the parameters of 
the poUing system are changed from cycle to cycle in a random manner. Since 
in the paper we are interesting in the characteristics related with busy periods 
of the system whose work starts by the arrival of a single customer to a station 
J e {1, ...,m} at moment 0, the law of arrival of customers to an idle system 
plays no role for subsequent arguments. 

The aim of the article is to study the distribution of the total size of the 
final product accumulated in the system during a busy period of the server. 
In particular, letting the final product (pi{n,j) be the service time of the j'-th 
customer served during the n-th visit of the server to station i, we provide 
conditions under which the tail distribution of the length of the busy period 
decays, as y — » oo, like const x y~'^ for some k> 0. 

Before we proceed to the rigorous statements of our results, consider two 
examples of branching type polling systems. 

Let 7+= {T} be the set of all probability distributions of nonnegative random 
variables, T™ := {(Ti, ...,Tm) : Ti e T^} be the set of all m-dimcnsional tuples 
of such distributions, A4e = {S} be the set of all w x m matrices £ — (£ij)™-^;^ 
with nonnegative elements, and Ai-y = {F} be the set of all m x (to+ 1) matrices 
^ — (7y)i^i j=o with nonnegative elements such that 

m 

^7ij=l, i=l,...,m. 
j=o 

Let P be a measure on the Borel cr— algebra of the space A4e x ^4-y x T^. 

ExEimple 1 (motivated by [34])- Consider a polling system with m stations 
and a single server performing cyclic service of the customers at the stations. 
Assume that initially there are no customers in the system and the server is 
located at parking place R. Assume that given the idle system the flow of cus- 
tomers arriving to station i is Poisson with, say, a deterministic rate Si . When 
the first customer appears in the system the server selects a random element 
(£:o,ro,To) e Me X X with 

£o= {eijmZ=i ' To = {lijm7=ij=o ' To = (Tio, T^o) 

and imm,ed,iately starts its zero service cycle (1 — > 2 — ^ • • • — > m) adopting 
the gated server policy with zero switchover times. Namely, the server serves 
all the customers that were queueing at a station when the server arrived and 
then instantly jumps to the next (in cyclic order) station. For the period while 
the server performs the batch of services at station i, new customers arrive to 
the system according to independent Poisson flows with intensities given by the 
vector (eii(O), £^2(0), eMn(O)) (some the components may be equal to zero) 
and the service times of customers are iid and distributed according to Tio{x) := 
P (7"i(0) < x). Each served customer either goes to station j e {1, . . . ,m} with 
probability 7y (0) or leaves the system with probability 7io(0) independently of 
other events. Besides, after the end of each service period of a customer the 
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customer contributes to the system its service time as the final product. It is 
assumed that given {Sq, Tq, Tq) the service times and the arrival process of new 
customers are independent. 

The subsequent routes n = 1,2,... have the same probabilistic structure spec- 
ified by the tuples 

with only difference that at the beginning of cycle n > 1 there is a possibility to 
have more than one customer in the system. 

Let us show that this system possesses a branching property. To this aim 
denote 

/•OO 

tiniX) ■■= / e-^^dTUx) 



the Laplace transform of the distribution Ti„(x) of the random variable Tj(n) 
representing the service time of a customer at station i during the n— th visit of 
the server. 

It is not difficult to check that the m.p.g.f. ^„(s;A) has the components (in 
our setting and the service time of customers as the final product) 



'(s;A) 



ii(n) Ji2in) 



>(")g--^Ti(") 



7io (n) + J2 '^J j n 



-Xx 



dTinix) 



i-i^e-^'^dTir^ix) 



7io(n) +^7ij{n)Sj \ tin I X + ^£ij(n)(l - Sj) 



(1) 



and 



/i«(s) := </.«(s;0) = 7io(n) + ^lijin)s, t„, ^ £,,(n)(l - Sj) . (2) 



i=i 



Example 2 (compare with [34]) ■ Consider the same polling system as earlier 
but assume now that at each station the server adopts the exhaustive server 
policy: it serves all the customers that were queueing at the staMon when the 
server arrived together with all subsequent arrivals up until the queue becomes 
empty and then instantly jumps to the next station. 

Let us show that this system possesses a branching property as well. 
Since during the n-th visit of the server to station i each customer at this 
station is served by the server a random number of times having shifted by 1 
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geometric distribution with parameter "fu^n), the Laplace transform Wi„(A) of 
the distribution of the random variable r]i{n), the total service time of a customer 
at station i during the n— th cycle, has the form 



(3) 



Let now ((Tii(n), . . . ,<Jim{n)) be a random vector distributed as the vector of 

the number of customers arriving to the stations of the polling system diiring 
the total service time of a customer at queue i which, after the end of its service 
at station i moves either to a station j ^ i ov leaves the system. Denote 



Vinis) ■■ = 



Similarly to (1) one can show that 



E 



CTii(n) _ _ _ „<7i,n{n) -XTi{n) 



yin{s)Wi„ j A + ^£y(n)(l - Sj) j . (4) 



Let now Im (A) be the Laplace transform of the distribution of a busy period 
'ili,tot{n') of the system generated by a single customer in an M/G/1 queue with 
arrival rate £i(n) and the service time distributed as r?i(n), and let {9i\{n), . . . , 9im{n)) 
be the total number of new customers arriving to the i-th stations of our polling 
system within the time- interval distributed as r]i^tot{n)- In this case lin{^) is a 
unique solution of the equation 

knW = UHn (A + Sii{n) (1 - /i„(A))) 

and (compare with [36]) the fmictions 



4^)(s;A) :=E 



e.i(n) e.2(n) 
Si S2 



,i = l,2,.. 



, m 



are (in our setting) unique solutions of the equations 



<i>l:\s;X) = yin{s)w,n A + ^£i,(n)(l - Sj) + £ii(n)(l - <^W(s;A)) J • 
In particular. 



/iW(s) = yin{s)win \ Y.eij{n){l - s^) + £i,(n)(l - ft«(s)) j • (5) 

Note that in the framework of the suggested approach one can consider 
also models with batch arrivals of customers as well. In this case one should 
replace, for instance, everywhere in Un{-) variables si, . . . , Sm by the respective 
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probability generating functions of the sizes of batches of customers arriving to 
stations j = 1, . . . , m. 

Since we imposed the branching property on the service discipUnes of the 
systems it is not a surprise for the reader that we eonduct our investigation 
by the methods of the theory of branching processes. This approach is not 
new. Polhng systems possessing the branching property in which the probabihty 
generating functions hn\s), n = 1, 2, . . . are nonrandom and the same for all 
cycles were considered in particular, in [24], [36], [42] and quite recently in [9], 
[12] and [13]. These models cover many classical service policies, including the 
exhaustive, gated, binomial-gated and their feedback modifications (see surveys 
[4] and [35] for definitions and more details). 

Polling systems with input parameters and service disciplines changing in a 
random manner and(or) depending on the states of systems are not studies yet 
in full generality. The analysis of such systems uses rather often the fiuid method 
([21]-[23]) or a method based on the construction of appropriate Lyapunov func- 
tions ([32] - [34]) (in the last case without reference to branching processes). 
The authors of [34] write that it is possible to generalize their models. However 
" ...This leads to a lot of complexities in the proofs." In the present paper we show 
that the reduction of the problems related to the branching type polling systems 
with final product and evolving in random environment (BTPSFPRE) to the 
respective problems for the multitype branching processes counted by random 
characteristics and evolving in random environment (MBPRCRE) gives the de- 
sired answers by a unique method and in a general situation. Since m.p.g.f.'s 
are selected for each cycle independently of the past, diff'erent service disciplines 
are allowed at different stations in our model. Moreover, the service disciplines 
at the stations may be changed at random from visit to visit. In particular, 
we allow for the mixture of exhaustive and gated service disciplines. Here we 
consider the models with zero switchover times. The study of BTPSFPRE with 
positive (and random) switchover times may be conducted by a similar method. 
This, however, requires more efforts and will be done elsewhere. 

The scheme of the remaining part of the paper looks as follows. Section 2 
is devoted to a detailed description of multitype branching processes counted 
by random characteristics and evolving in random environment. In Section 
3 wc recall some known results for ordinary multitype branching processes in 
random environment (MBPRE) and formulate the main results of our paper 
describing the asymptotic behavior of the total size of the final product for 
subcritieal MBPRCRE. Section 4 recalls an important statement related with 
the asymptotic properties of the tail distribution of infinite sums of products of 
random matrices. In Section 5 we deduce some estimates for the moments of the 
population sizes of subcritieal MBPRE's. The proofs of the main results of our 
paper are collected in Section 7. And, finally, we demonstrate in Section 8 how 
one can use the obtained results for MBPRCRE to make conclusions about the 
probabilistic properties of various characteristics of the branching type polling 
systems. 
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2 Branching processes in random environment 
counted by random characteristics 

As we claimed in the previous section, the main goal of the present paper is 
to analyze properties of busy periods and some other cliaracteristics of the 
brandling type polling systems by means of MBPRCRE. However, to start 
such analysis we need to pass through a relatively long way of notation and 
statements which is practically always the case when one consider multitype 
branching processes. 

Let {^1, ■ ■ ■ ,^m',^) be a (m + l)-dimensional vector where the components 
^1, . . . , are integer- valued nonnegative random variables and is a nonneg- 
ative random variable and let 



F{s;X) = E 



*1 •'2 



* 5™ 6 



, s = (si, . . . , Sm) e [0, 1]™ , A > 0, 



be the respective m.p.g.f. Denote := {F(s;A), s € [0, 1]™} the set of all such 
m.p.g.f.'s and let 

:=J'xxJ'xX---xJ'x= {f(s;A) = (f^'\s;X), . . . ,F(™)(s;A)) } 

be the m-times direct product of J^\. Let, further, 

.Fo:={/(s) = i^(s;0),se [0,1]"} 
be the set of all ordinary probability generating functions (p.g.f.'s) 

and 

^™ := ^0 X • • • X .Fo = {f(s) = (/(i)(s), . . . , /(-)(s)) } 

be the set of all m-dimcnsional (vector- valued) p.g.f.'s. Assume that a probabil- 
ity measure P is specified on the natural a-algebra A generated by the subsets 
of JPj^. Let 

Fo(s;A), Fi(s;A), . . . , Ffe(s;A), . . . with F„(s;A) := (f^^) (s;A), . . . , F(r\s;X)) 

be a sequence of vector- valued m.p.g.f.'s selected from .7^™ in an iid manner in 
accordance with measure P . The sequence {F„(s;A), n > 0} is called a random 
environment. With m.p.g.f. Fn^ (s; A) we associate a random vector of offsprings 
^i(n) =: {^ii{n),^i2{n), . . . ,£,im{n)) and a random variable ipi{n) such that 

F«(s;A) = E ^("^ • ■ • .se.™(n)g-Av,,(n)" 

Now we may give an informal description of the MBPRCRE 
R(n) = (Z(n);$(n)),n = 0,l,..., 
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which may be treated as the process describing the evolution of a m— type 
population of particles with accumulation of a final product. 

The starting conditions of the process are: a vector (may be random) of 
particles Z(0) = (Zi(0), . . . , Z„i(0)) where Zi{0) denotes the number of particles 
of type i G {1, . . . , m} in the process at moment 0, and an amount $(0) (may 
be random) of a final product. All the particles have the unit life length and 
just before the death produce children and final prodiicts independently of each 
other. For instance, a particle, say, of type i produces particles of different 
types and adds some amount of the final product to the existing amount of 
the final product in accordance with m.p.g.f. _Fq*''(s;A). The newborn particles 
constitute the first generation of the MBPRCRE, have the unit life-length and 
dying produce, independently of each other, offsprings and final products in 
accordance with their types and subject to the m.p.g.f. F^^^ (s;A), i = 1, 2, . . . , m 
and so on. 

A rigorous definition of the process we are interesting in looks as follows. 

Definition 3 A (m + I) -dimensional Galton-Watson branching process 

R^(n) = R(n) := (Z(n); $(n)) = (Zi(n), . . . , Z„(n); $(n)) , n £ No, 

counted by a random characteristics (p in a fixed (but picked at random) en- 
vironment {F„(s;A),n >0} is a time-inhomogeneous Markov process with the 
state space 

xKo:= {z = {zi,...,Zm;w), ZiGNo;wG [0,oo)} 

defined as 

R(0) = (Z(0);a>(0)) = (z;(^o), 
E [s2(n+i)e-A*("+i) I Fo, . . . , F„; R(0), . . . , R(n)] = e-^*(") (F„(s;A))^(") . 

(6) 

Note that the initial value R(0) may be random, and, for the reason of 
applications to queueing systems we do not exclude the case z = 0. 
In what follows, to simpliiy notation, we write 



Thus, R(0) = (z;93o), 



gZ(n+l)g-A*(n+l) 



E 



s2("+i)e-^*("+i)|Fo,...,F„ 



Ep 



gZ(n+l)g-A$(„+l)| R(0),,...,R(n) 



= e 



(F„(s;A))^(") 



or 



Z.(n) 



R(n+ 1) 
$(n+ 1) 



(0;$(n)) + ^ J2 {Un\k);'Pi{n;k)), 

i=l k=l 
n m Zi{l) 



(7) 



;=o i=i fe=i 
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where the random vector (^i(n; fc); (pi{n: k)) represents the offspring vector and 
the size of the final product of the k—th particle of type i of the n— th generation 
of the process. Given the environment and n = 0, 1, . . . , and i G {1,. . . , m} the 
vectors 

(Ci(n; k); (pi{n; fc)) , fc = 1, 2, . . . , Zi{n) 

are independent and identically distributed: (^i(n; k);ipi{n; fc)) = (^i(n); tpi{n)). 

Observe that if $(0) = and (pi{n;k) = 1 then $(n) is the total number 
of particles born in the process within generations 0, 1, n — 1; if (pi{n;k) = 

J 1^^-^ (n; fc) > t| for some positive integer t (here and in what follows 

/ {A} means the indicator of the event A) and $(0) = 0, then is the total 
number of particles of all types in generations 0, 1, n — 1 each of which had 
at least t children, and so on. 

Letting A = we arrive to the definition of the ordinary multitype branching 
process in random environment (MBPRE) which we call the underlying MBPRE 
for the initial MBPRCRE. 

Definition 4 A m-type Galton- Watson process 

Z{n) = {Zi{n),...,Zrn{n)),neNo 

in a fixed (but selected at random) environment {f„(s), n >0} is a time-inhomogeneous 
Markov chain with the state space 

N™:={z=(^i,...,2„),z, gNo} 

defined as 

Z(0) = z, E [s^("+i) I fo, . . . , f„; Z(0), . . . , Z(n)] = (f„ (s))^^") . (8) 
To simplify notation we write 



Ef 



,Z(«+1) 



E 



s^("+i)|fo,...,f„ 



It follows from (8) that Z(0) = z and 



m Zi(n) 

z(n+i) ■.= Y.Y. ^'('^;^) 
i=i fe=i 

where ^i{n;k) = ^i(n),fc = 1, 2, . . . , Zj(n), and, given the environment and 
n = 0, 1, . . . , and i G {!,..., m} the mentioned random vectors are independent. 
Let 



A„ = (aiiH)™. 1 := 



|s=l,A=0 



dsj 



(9) 
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be the mean matrix of the vector-valued p.g.f. fn and 

C„ := (Ep^i(n),...,EpvPi(n))' = (^^^^^ |s=i,A=o , ^^^^^^ |s=i,A=o 

(10) ' 

By our assumptions the pairs {An, C„) , n = 0, 1, ... are iid: {An, C„) = {A, C). 
Suppose that 

Elog+ ||A|| < cx). (11) 
It is known (see, for instance, [31]) that given condition (11) the Umit 

lim -log\\An-iAn-2---Ao\\=:a (12) 

n— *oo n 

exists with probabihty 1 and, moreover, 

Um -E log \\An-iAn-2 ■■■Ao\\=a. (13) 

n— »oo n 

In what follows we call a MBPRE subcritical if a < and supercritical if a > 0. 

Recall that the single- type BPRE with iid offspring p.g.f. 's were introduced 
by Smith and Wilkinson in [37] and, in a more general setting, in [14] , [15], 
[16] and have been investigated by many authors (see survey [39] for a list of 
references up to 1985 and [1], [2], [3], [7], [10], [11], [20], [25], [40] and [41] for 
some more recent results). MBPRE were analyzed, in particular, in [14], [30] 
and [38]. 

Ordinary single- type Galton- Watson branching processes counted by random 
characteristics whore investigated by Scvastyanov [8] ( for integer-valued ip{n: k)) 
and by Grishechkin [6] (for the general (f{n;k)). Grishechkin used the Galton- 
Watson and continuous time Markov branching processes counted by random 
characteristics (with or without immigration) to study queueing systems with 
processor sharing discipline [5]. 



3 Limit theorems for MBPRCRE 



Introduce the notation 

n-l 

Tll,n := n 1 < ^ < 

i=l 

with the agreement that Iln,n '■= E is the unit m x m matrix. For vectors 
u = {ui, . . .,Um),'v ={vi, . . . ,Um)'eR™ denote 

m 

(u> v) := ^ UiVi 
fe=i 
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their inner product. 

For &mxm matrix A = {aij)^j_^ and a m— dimensional vector u =(ui, . . . , Wm) 



introduce the norms 



||A||:=5^K|, ||u||:=^ 

i,j=l i=l 



and 



11^112 ■■ = 



\ i=l 



Now we formulate an important statement concerning properties of MBPRE. 
Let 

Qi (f) := hm Pf (||Z(n)|| = 0|Z(0) = 

n— *oo 

be the extinction probability of a MBPRE initiated at time by a single indi- 
vidual of type i and 

q(f) := (91 if),..., qm{f)). 

Theorem 5 ([S8J) If the mean matrices of a MBPRE meet condition (11) and 
there exists a positive integer L such that 



{ min {AL-lAL-2■■■AQ)^.>{y\=l 

yl<i,j<m ■' J 



and 1 < I <m such that 

E |log (1 - Pf (Zi(L) = OlZ(O) = e,))| < 

then, for a specified by (13) 

1) a<0 implies P (q(f) = 1) = 1; 

2) a>0 implies P (q (f ) < 1) = 1 and 

Pf ( lim log ||Z(n)l| = a|Z(0) = e,) = 1 - (f) 
with probability 1 for 1 < i < m. 



(14) 



Let T be the extinction moment of a MBPRE which starts by (may be 
random) vector of the mmibcr of particles Z(0) with E ||Z(0)|| < oo. Clearly, 

P(r>n) = P(||Z(n)|| > 1) <E||Z(ri)|| =E||Z(0)^oAi---^„_i|| 
< E||Z(0)||EPoAi---A„_i||. 

Hence we see that if a < then for any a, G (0, —a) there exists a constant 
if* = ii'*(a*) G (0, oo) such that for each n = 0, 1, 2, . . . 



P (r > n) < isT.e-"*". 



(15) 
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Let 

oo m Zi (n) 



$ := lim $(n) = $(0) + V V V ^i(n; k) 

n— »oo -^^ — ' ^ — ' ^ — ' 



n=0 i=l fc=l 



be the total size of the final product produced by the particles of the MBPRCRE 
up to the extinction moment (if any). It is easy to see that if the underlying 
MBPRE is supercritical, satisfies conditions of Theorem 5 and, in addition, the 
final product of the MBPRCRE meets the condition 



P min Ef(^i(n)>0 >0 (16) 
then, by (14) and the law of large numbers, for each i = 1, . . . ,m 

P ($ = (X)|Z(0) = ei) > P lim Jnf^ ^ ^ ^j(n; k) = cx)|Z(0) = eA > 0. 

(17) 

Note, finally, that if $(0) = and ^j{n] k) = \ then 

oo 

$:= lim $(n) = y ||Z(n)|| 

is the total number of individuals ever existed in the MBPRCRE or, what is 
the same, in the underlying MBPRE. 
For a given a; > set 



s{x):= lim (E||A„_i-..Aory''"= lim (E ||no,„|r)'/" (18) 

n— >oo n— »oo 

and let 



s'(0)= lim -ElogP„_i---Ao|| = lim -Elog ||no,„|| (19) 

n — ^oo Ji n — >QO n 

be the top Lyapunov exponent for this sequence of matrices. 

Denote I> := {x > : E ||Ao||^ < oo}. It is known that the limits in (18) and 
(19) exist and, moreover, s{x) is a log-convex continuous function in D (see, for 
instance, [31]). Put 

K :=inf {a; > : s(a;) > 1} (20) 

and K = 00 if s{x) < 1 for all x > 0. Observe that s(0) = 1 and, therefore, 
K = if s'(0) > and K e (0, oo] if s'(0) < 0. 

In the last case (which will be our main concern) the series YL^=o ^ ||no,nir 
converges if < a; < «; and diverges if a; > k. 

Introduce the set 

U+ = {a= (ui, . . .,Um) e : Mi > 0, 1 < i < m, ||u||2 = 1} 



12 



and associate with the tuple {An, C„), n = 0, 1, 2, ... of iid pairs the series 

oc oo 

El := ^'^AiAi+i . . . A/c_iCfe = ^Ilj^feC/c, I = 0,1, . . .; S =: Eq. 

k=l k=l 

Our main results are established under the following hypothesis. 
Condition T. There exist positive constants k and Kq and a continuous 
strictly positive function /(u) on U+ such that for all u Gt/+ 

lim y«P((u,S> >2/)=i^o/(u). 

In Section 4 we list sufficient conditions on the distributions of the pairs 
{An, Cn) which provide the validity of Condition T. These conditions are ex- 
tracted from paper [28] where the behavior of the tail distribution of sums and 
products of random matrices were investigated. 

The following theorem is the main result of the article. 

Theorem 6 Let a MBPRCRE satisfy the following hypotheses: 

1 ) the underlying MBPRE is subcritical and meet conditions of Theorem 5; 

2) for K specified by (20) the following assumptions fulfill: 
if K> 1 then 



max E 

l<i<m 



< 00 and E 



^{ipi{n) -'E,FiPi{n)) 



< 00, (21) 



if K <1 then 

I m 

^max E I VavF^ij j < oo and E I VarFifii{n) J < oo; 



(22) 



4) there exists 6 > such that < Etp'^^^{n) < oo, i = 1, . . . ,m. 
If, in addition, the mean matrix (9) and the vector (10) are such that Con- 
dition T is valid, then, as y ^ oo 

P($>y)~C2/-^ Ce(0,oo). 

An evident corollary of Theorem 6 is the following statement. 

Theorem 7 // the conditions of Theorem 6 are valid then 

E^'^ < 00 

if and only if x < k. 
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4 Auxiliary results 



The proof of Theorem 6 is heavily based on Condition T whose vahdity is not 
easy to check. We list here a set of assumptions given in [28] which imply 
Condition T. 

Let A (A) be the spectral radius of the matrix A. The following statement 
is a refinement of a Kesten theorem from [28] . 

Theorem 8 (see [18]) Let {An,n > 0} &e a sequence of iid matrices generated 
by a measure with support concentrated on nonnegative matrices and A = 
i^ij)Tj=i ^ Assume that the following conditions are valid: 

1) there exists e > such that E ||^||^ < oo; 

2) A has no zero rows a.s.; 

3) the group generated by 

{log A(a„ • • • ao) : a„ • • • ao > for some n and Ui G supp{]Pa)} 

is dense in R; 

4) there exists Kq > for which 









E 


min 

l<i<m 











and 

EllAf log+||A|| <oo. 
Then there exists a k € {0, kq] such that 

s'{k)= lira ilogE||A„_i---^oir = 0. 

n^oo n 

If, in addition, the tuple of m- dimensional vectors {C„,n > 0} is such that 
the pairs (A„, C„), n = 0, 1, . . . are iid: {An, C„) = {A, C) and such that 

P(C = 0)<L P(C>0) = 1, E||Cf<oo, 

then there exist a constant Kq G (0, oo) and a continuous strictly positive func- 
tion l{vL) on U+ such that 

hm y-P ((u,S) >y)= KoKu), u £[/+. 
y—KX 

The next lemma will be of importance for subsequent arguments. 

Lemma 9 ([27], Theorem 1.5.1) If Xi, i — 1,2, ... is a sequence of iid random 
variables such that E \Xi\^ < oo and EXi = ifp > 1, and N is a stopping time 
for the sequence Sn = Xi + . . . + X„, then there exists a constant Rp G (0, oo) 
such that 

E\SNf < RpE\XifENP/^''\ 
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5 Properties of the underlying MBPRE 



In this section we assume that the conditions of Theorem 6 are vahd. This 
means, in particular, that we deal with subcritical MBPRE. Let us agree to 
denote by K,Kx,x G (0, oo) positive constants which may be different from 
formula to formula. 

First wc evaluate the expectation of the random variable ||Z(n)||^ , < x < 
K, from above. 

Lemma 10 //Z(0) = z and k > 1 then for each x G [1, k) there exist px G (0, 1) 
and Kx < oo such that 

E||Z(n)r <^r.p:j||zr (23) 

for all n = 1,2, 

Proof. Clearly, for any nonrandom vector b € M"* 

E(Z(n),b) = eIJ2 E ei(n-l;fc),b\ 

\i=l k=l I 

= E (Z(n - 1), An-ih) = E (Z(n - j), n„_,>b) 
= . . . = E (z,no,„b) , < j < n. 

Thus, if ||b|| < K < (X> then for any S > with s(l) (1 + 15) < 1 there exists a 
constant Ki G (0, oo) such that for all n = 1, 2, . . . and < j <n 

|E(Z(n-i),n„_,-„b)| < ||z||E||no,„||E||b|| <i^i(s(l)(l + 5))"||z||. (24) 

Now we apply arguments similar to those used in [19]. It is easy to check 
that for any y,w >0 and any e G (0, 1) 

t/- < {1 + S)W- + Cx,e\y-W\^ , 

where Cx,e '■= — (1 + s)"^^^) • Hence we have 

E||Z(n)|r = E||Z(n-l)A„_i+Z(n)-Z(n-l)A„_if 

< (1 + £)E \\Z{n - l)A„_i f + Cx,eE ||Z(n) - Z(n - l)A„_i f . 

(25) 

Recalling the definition aij{n) = EF^ij(n), set 

m 

Piin):=J2^^ij{n)-aij{n)) (26) 



and let 



Mx{n;i) := Ep 1^(71)1"^ , Mx := max E |A(n)r • 

l<i<m 
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By Lemma 9 with p = a; > 1 we conclude 
Ef ||Z(n)-Z(n-l)A„_if 

rn Z^{n-1) rn 

=^^12 12 12 l^iJ in-l;k)- aij (n - 1)] 

i=l fc=l j=l 



Zi(n-l) m 

^^ij {n-l;k)- aij {n - 1)] 

/c=l j=l 



< ^ Ef 

i=l /c=l j = l 

m 

< ii^m^ ^ M^(n - 1, i)EFZf /'""^n - 1). 



(27) 



These estimates and (25) give 

E||Z(n)|r < (l + £)E||Z(n-l)A„_ir 

m 

+c^,,R^m^M^ J2 ^Zf^''\n - 1). 



(28) 



Now we are ready to demonstrate (23) . First we estabhsh (23) for all integer 
X e [1, k). For a; = 1 we have proved (23) by (24) with b = 1. Now we use 
induction for x >2. Observing that x/2 V 1 < a; — 1 in this case, we see by (28) 
and the estimate 

m m 

^EZf/^^^n - 1) < E^ Zr^(n - 1) < E ||Z(n - l)\r' 

i=l i=l 

that 

E||Z(n)|r < (l+£)E||Z(n-l)^,_ir 

+c^,eta"Af,E||Z(n- l)|r"' < (l+e)2E||Z(n-2)^„_2^„_i 
+ {l + £)cx,eR.m-M^E WAn-if E \\Z{n - 2)\r' 
+c^,,R^m''M^E \\Z{n - l)f < • • • < (1 + £)"E ||zno,„||^ 

n-l 

+c,,,R,m^M, ^(1 + ey-E ||n„_,- „|r E ||Z(n - j - l)f-^ . 

3=0 

Since x G [1, k), for any S > there exists a constant such that 

Eiino,„ir <L. {s{x){i+s)r 

for all n = 0, 1, 2, .... By induction hypothesis there exist constants K^-i and 
Px-i G 1) such that 



E ||Z(n - j)\r' < K^-^plZi < i^x-ipH 11^ 



"-3 



16 



for all j = 0,1, . . . ,n. Thus, 

E ||Z(n)r < (1 + erz, {s{x){l + 5))" Hzf 

n-l 

+c,^,R,m-M,L,K,_^ Hzf ^(1 + ey {s{x){l + 5)y p^zf 

=L,s^ix){i+6r{i+5r\\zr 

n j—1 

Now selecting S and e in such a way that 

S{x)il+S) (1 + 5) G (Px-1,1) 

(the last is always possible since p^-i < s{x) < 1) we get (23) . 

To treat the case of noninteger x € [1,k) observe that x* = [k] G [1,^) is 
an integer for which (23) is vahd. Thus, it remains to demonstrate (23) for 
X = X* + J < K, where 7 G (0, 1). Since x* > x/2 V 1 one can use the same 
arguments as earher with x* for x — 1. 

The lemma is proved. 

Let r be an integer and 

C = C('') •= rninjn > : ||Z(n)|| > r} 
with the natural agreement that = 00 if max„ ||Z(n)|| < r. 
Lemma 11 Under the conditions of Theorem 6 for any fixed r > 1 
E[||Z(C(r))|rJ{C(r)<oo}]<oo. 
Proof. Similar to (28) we have for a; > 1 

Ef [||Z(n)|r |Z(0), . . . , Z(n - 1)] < (1 + e) ||Z(n - l)A„_ir 

m 

+c^,eRxm'' J2 Mx{n - 1, i)Zf /'""^n - 1) =: 4-^(71 - 1), (31) 

while Jensen's inequality yields for a; < 1 

EF[||Z(n)r|Z(0),...,Z(n-l)] < (Ep [||Z(n)|| |Z(0), . . . , Z(n - 1)])^ 

= ||Z(n-l)A„_i||" < *^(n-l). (32) 

Clearly, 

*«(n-l)/{n<C} < Q„-i(r):=(l + £)r'^m„_i|r 

1=1 
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Using these estimates we have on the event {C < 00} : 

l|z(C)r = *«(C-i)^^^ < Qc-iM "^^^^"'^ 

By (31) and (32) we now see that 

E[||Z(C)ir7{C<oo}]<5^E 



(33) 



C<n<T 



n>l 



< ^ E [Qn{r)I {r > n}] = ^ P (r > n) E [0„(r)] 

n>l n>l 

= [(1 + e)r^E \\Af + c,,,m^+'r^/^''^M,'\ ^ P (r > n) < 00, 

n>l 

since for m := min (k/2, 1) 

X^P(r>n) = ^ P (||Z(n - 1)11 >l) = ^E[PF(||Z(n- 1)11 >1)] 

n>l n>l 

< ^E[EF||Z(n-l)f^] < ^E||EFZ(n-l)|r^ 

n>l n>l 

< ^iizr^E||Ao---A„_ir^ <oo. 



n>l 



n>l 



The lemma is proved. 

Let S„ (n = 1, 2, . . .) be the cr-algebra generated by the tuple 



and let 



Fo(s, A), Fi(s, A), . . . , F„_i (s, A), Z(0), . . . , Z(n) 



C„ := (EF(/?i(n),...,EF¥'m(n))', n = 0,1,2,... 



Recall that the pairs {A^, C„) , n = 0, 1, ... , are iid according to the defi- 
nition of our MBPRCRE and, in particular, {An, C„) is independent on B„ . 
Set 

00 

5(C):=$^(Z(n),C„). (34) 

The next lemma shows that for large r the random variable S{Q = S{C,(r)) 
is, in a sense, close to the conditional expectation E [5(^)1 S^] = (Z(C), S^) . 

Lemma 12 Under the conditions of Theorem 6 for any e > there exists 
r = r(e) such that for all y >yo 

P(|5(C) - (Z(C),Sc)| > £y;C < cx)) < ^E [l|Z(C)r / {C < 00}] . (35) 
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Proof. Evidently, for n > C + 1 

n 

(z(n)-z(c)nc,„,c„) = (z(Onz,„-z(;-i)n,_i,„,c„) 



;=c+i 



which impUes 

|5(C)-(z(C),Sc)| 



= J2 (z(0-zG-i)^;-i,nz,„c„) 

i=C+i 



^ (Z(n)-Z(C)nc,„,C„) 

n=C+l 

CO 

< Yl l(Z(n)-Z(C)nc,„,C„> 



n=C+l 



E 



J2 (z(0-z(/-i)A;_i,ni,„c„) 



;=c+i 



< Yl \\m-^{i-mi-i\ 



Y ||z(0-z(z-i)^,_i||||s,||. 



(36) 



Hence, on account of J2p=ij ^ — ""^/S < 2, we may apply on the event 
{C < 00} the arguments used in [29], Lemma 3, to conclude that 

Pp(|5(C)-(Z(C),S^)|>ey|Bc) 

(00 00 ^ \ 

Y ||z(0-z(/-i)A,_i||||s,||>6rty Y TTTT^I^c 

< Y Pf(||Z(0-Z(Z-1)A,_i||||S,||> 



2(z-cr 



(37) 



Since Z{1) — Z{1 — l)Ai-i and are independent random vectors on the event 
C < Z < 00, we get 



Pf \\Z{l)-Z{l-l)Ai_,\\\\Ei\\ > 



ey 



2{l-Cr 



= Pf (||Z(0 - Z{1 - l)Ai_4 ^dt\Bc)V {^\E.\\ 



> 



2t{i-cY 
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According to Condition T there exists a constant K G (0, oo) such that for all 



Pf \\Z{l)-Z{l-l)Ai_i\\\\'Ei\\ > 



ey 



< 



2(1-0 

/■oo ly-.K 

/ Pf (||Z(/) - Z(Z - l)A-i|| G dt |Bf ) -— {I - C) 

Jo ^ y 



<^A^- 0'" Ef [||z(0 - z{i - I) A., r \Bc]. 
^ y 

Now we consider the cases k < 1 and k > 1 separately. 
For the first case we use the estimate 

Ef [||Z(0 - Z{1 - l)A-iir |Sc] < (Ef [||Z(0 - Z{1 - l)Ai_if \ Bf ]) 

Further we have 
Ef 



(38) 



k/2 



(39) 



iiz(/)-z(/-i)Ai_iinsc 



E E Y.\-^ij{l-'^\k)-ai^{l-l)]] \Bc 
i=i fc=i j=i 



= ^EF/3^(i-l)EF [Z,(Z-l)|Sc] 

m m 

= E^F/?- (« - 1) (z(C)nc,;-i), < ||z(c)nc,;-i||EEF/?f(; - i). 



i=l 



i=l 



(40) 



Thus, for «; < 1 



(ef [||z(0 - z{i - i)Ai.,f \Bc]y^' < wzma-if/' ( E^F/?! (« - 1) 



k/2 



(41) 



This, in view of the inequality 



s{k/2) = lim (Ellno.nf/') 



1/n 



< 1, 
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the first part of condition (22), and relations (37)-(41) leads to the estimate 
P(|5(C)-(Z(C),Sf)|>£2/;C<oo) 

K 



< 



K 



-E 



^ a-c)'niz(c)ir/'iinc,i_ir/'/{c<^} 



i=f+i 



-E 



|z(c)ir/'/{C<oo}X;«"'E||no,,-i| 



k/2 



1=1 



< ;;^^;-3^E[||Z(C)ir/{C<oo}]<i,E[||Z(C)ir7{C<^^^ (42) 

for all r > ro(e), proving the lemma for «; < 1. 

For the case k > 1 we use Lemma 9 to conclude that for any I > ( 

EF[||Z(Z)-Z(/-l)A-ir|Bd 

m 

< R,m''Y.M,{n;i)Ep [\\Z,{1 - 1)11"/"^^ \Bc 

i=l 

m 

< i?«m«Ep [||Z(/ - 1)11"/^^^ I Sc] Y,M.in;i). 

i=l 

By Lemma 10 there exist constants Pre/2vi € (0, 1) and K^/2\/i < oo such that 
for alU > C 



E. 



|Z(Z-l)r/2^i|Bc] <i^«/2vipL7^vi'l|Z(C)ll 
This yields the estimates 

P(|^(C)-(Z(0,Sc)|>ey;C<oo) 



i/2Vl 



< 



< 



-E 



/{C<oo} 



-E 



E a - C)'" Ep [ II Z(/ - 1) f/^^^ \bAi{C<<x>} 



i=C+i 



^KyK 



J2 (i- 0"= iiz(or/'^' i^./2vipL72vi^{c < oo} 

i=C+i 



E 



z(c)r/^"^/{c<oo}E^^>L72 



/2V1 



< ,.,C"-!/2vi E[iiz(c)rj{c < 00}] < ^E[||z(c)r /{c < oo}] 

(the last is valid by selecting r sufficiently large) which justifies the statement 
of the lemma for k > 1. 
The lemma is proved. 



21 



6 The accumulated amount of the final product 

In this section we deduce some estimates related with the total size of the final 
product accumulated in a subcritical MBPRCRE during its evolution. 
Let 

m Zi (n) 

i=i k=i 

be the total size of the final product produced by the individuals of the n-th 
generation of a MBPRCRE and 

oo 

^N) := Mn), * = ^0- 

ra=JV 

Lemma 13 Under the conditions of Theorem 6 for any e > there exists 
r = r(e) such that for all y > yo 

P (|i(C) - 5(01 > £y;C < oo) < ^E[||Z(C)r /{C < oo}] . (43) 
Proof. We have 

m Zi(n) 

A(n)-(Z(n),C„)=^ ^ (y)i(n; A;) - EF^i(n)) . 

i=i k=i 

Now we consider separately the cases k < 1 and k > 1. 
For K > 1 we use Lemmas 9 and 10 to get for n > C : 



Pf |A(n)-(Z(n),C„)|> 



2(n-C)2 



4-(n-C)^ 

- 1 V5 



^ ^ (!/?,(n;A:) -EF<y£'i(n)) 

j=i fe=i 



< (4m)'=i?«^^y^^EFki(n)-EF^i(n)rEF ||Z(n) 



2k ™ 



{ey) 



i=l 



ik/2V1 



< (4m)-i?,i^,/2vi^^^^^^ J^Ef |¥'i(n) - EF<^,(n)r ||Z(C)|| 



k/2V1 n-C 

^k/2V1' 
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Hence, in view of condition (21) 

P(|$(C)-5(C)| >e2/;C<oo) 



< E 



< 



L"=c 



^Pf (^IA(n)- (Z(n),C„)| > ^^-f^ |Sc j /(C < ^) 



E 



E 



oo 

5(C)ir/^"^/(C<oo)]^n^X 



/2V1 



< 



COT) r' 

(,,)«,.-(./.vi) E[iiz(c)r/(c < oo)] < -E[||z(c)ir/(c < oo) 



for all r > r{e). 

To analyze the case k < 1 we apply for n > C the inequality 



k/2 



Ef [|A(n)- (Z(n), C„)r l^d < (Ep [|A(n)- (Z(n), C„)|' | ]) 
Further, we have 

Ef [|A(n)- (Z(n),C„)|' 1^^] < ^Ef/S.'HEf [Z,(n) |S. 

m m 

= m2^EF/3^(n) (Z(C)nc,„), < ||Z(C)nc,„|| ^ EfA^^. 
Thus, for K < 1 

(Ef [|A(n)-(Z(n),C„)|^|ec ) ' < m« ||Z(C)nc,„+i f/^ ( ^ EfA^H 



i=l 



«/2 



k/2 



This combined with the assumption (22) shows that 
P(|$(C)-5(C)| >£y;C<oo) 



< E 



< 



^ Pf (^I A(n)- (Z(n), C„)| > ^^^^^^ I C, Z(0), . . . , Z(0 j / (C < oo) 



E 



J2 (/-C)"^l|z(C)ir/'||nc,;ir^/(c<oo) 



-E 



z(c)r/'E^"'E||no,Hr/'/(c<oo) 



< ;;^^E[||z(c)ir/{c < oo}] < -E[||z(or /{c < -}] 
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for all r > r{e). 

The lemma is proved. 

Up to now we have assumed that the initial number of particles and the 
initial volume of the final product are nonrandom. Lemmas 14 and 15 are free 

of this restriction. 

Lemma 14 Let a MBPRGRE be subcritical, E||Z(0)|| < oo and there exist 
6 > such that 

E$«+'5(0) < 00, max Eip'^+^ < oo. 

l<i<m 

Then for any e > there exists r = r(e) such that for all y >yo 

P($(C-l)>ey:C<oo)<-^. (44) 

Proof. Let, as earlier, r be the extinction moment of the underlying 
MBPRE and let c be a constant such that for all y > yo{c) 



t>clny ^ 

Such a constant, clearly, exists in view of (15). Recalling (7) put 

n m r 

$.(n) = $(0) + EEE^*(^;^)- 
;=o i=l k=l 

Now we have for y > yo{c) 

oo 

P ($(C - 1) > £2/; C < oo) = ^ P ($(n - 1) > C = n) 

oo oo 

< ^ P {^r{n - 1) > ey; C = < ^ P {^r{n - I) > ey\T > n) 

n=l n— 1 

n m r ^ x 



+ E P (r > n) < -— ^E$«+^(0) 

n>c\ny \ tJ) 

\^») l<n<clny 1=0 i=l k=l ^ 

< K ^ + — < 



Lemma 15 Under conditions of Lemma 14 for any r there exists yo = yo{r) 
such that for all y > yo 

P($>y;C = oo)<-^. (45) 
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Proof. Letting c be the same as in the previous lemma we have for y > yo{c) 



P($>y;C = oo) <P($^(r) >t/) <P($^([cln?/]) >y) + P(T>clny) 

'z\ny] rn r 



< P(<I>(0)>|)+ ^ ^^PU',(^;fc)> 

1=0 i=l fe=l ^ 



y 



2mr [c In y] 



+P (r > c\ny) < K 



In 



K+5+1 , 



i/K+l5 



i,K + l5 



as desired. 



7 Proof of Theorem 6 

Now we are ready to prove the main result of the paper, Theorem 6. First ob- 
serve that by the equivalence of the norms ||-|| and ||-||2 and estimates (35), (43), 
(44) and (45), for any e € (0,1/3) one can find r = r{e) such that for all 
y > yo{r,e) 

P($>y)<P($>y;C = oo) 

+P ((Z(C), Sc) > y(l - 3£); C < oo) + P (|$(C) - 5(C)| > ey: C < oo) 
+P (|5(C) - (Z(C), Sc)| > £y; C < ^) + P (1>(C - 1) > ey; C < ^) 

< P ((Z(C),Sc) > 2/(1 - 3£);C < 00) + [||Z(C)||^ / {C < oo}] + 

(46) 

Let K{1) :— infueu+ '(u) > where 1{vl) is the function involved in Condition T. 
By this condition and the independency of Z(^) and we conclude 



lim sup tz-^P ((Z(C), S^) > 2/(1 - 3e); C < oo) 



< lim sup 2/"/ P(Z(C) Gdu;C<(X))p((^,S 

S/^oo J||u||2=r V\l|U|l2 

= Ko{l-3e)-'^ r P(Z(C)erfu;C<oo)||u||^z(^ 

J\\u\L=r \\m\: 



> 



(l-3£)y \ 
l|u||2 J 



= Ko(l - 3e)-''E 



|Z(C)ll2M^^)^{C<oo} 



< 00. 



Since y'^P ($ > y) does not depend on r and s, the previous estimate and (46) 
yield 

lim sup y^P ($ > 2/) < 00 (47) 

y—>oo 



and, moreover, 

lim sup 2/'*P > y) < Kq lim E 



m)\\2i(j^^]Hc<oo} 



■ (48) 
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To get a similar estimate from below we use for e > the inequality 

P > y) > P (<& > y; C < 00) > P ((Z(C), Sc) > y(l + Se); C < oo) 
-P (|$(C) - S{0\ > ey; C < oo) - P (|5(C) - (Z(C), S^)! > sy;C < oo) 
-P ($ (C - 1) > ey; C < oo) . 

Now we select r as large to meet estimates (35), (43) and (44). This gives for 
sufficiently large y > r the inequality 

P ($ > y) > P ((Z(C), Sc) > y{l + 3s); C < ^)-^E [||Z(C)||2^ / {C < oo}]-^^. 

Letting y ^ oo we obtain 

lim inf 2/'^P((Z(C),S^) >2/(l + 3e)) 

y—^oo 



lim inf j/" 

y—^oo 



P(Z(C) e q'u;C < oo)P 
Z(C) 



" s\ > (^ + ^^)^ 



u 



u 



l|Z(C)ll2^ 



7{C<oo} 



1Z(C)II2 

We know that for sufficiently small e > and an appropriate r 

m 



iro(l+3£)-'^E 



|Z(C)II: 



J{C<oo} 



1|Z(C)II^ 

> {KoK{l){l + 3£)-« - 2e) E [||Z(C)||^ /{C < oo}] > 
which implies 



-2£E[||Z(C)||2^/{C<oo}] 



lim inf y^P ($ > y) > 



leading in turn to 



lim inf y^V ($ > y) > lim E 

y^oo r — ^■oo 



l|Z(0ll2' 



z(C) 



|Z(C)II: 



/{C<^} 



(49) 



This combined with (48) gives 
lim y'^P ($ > y) = Kq lim E 

y^oo r— >-oo 



l|Z(C)ll2i 



z(C) 
llz(C)lh 



/{C<oo} 



e (o,oo). 



The theorem is proved. 

Theorem 16 Let conditions of Theorem 6 he valid for a subcritical MBPRCRE 
starting at moment Q by a random tuple Z(0) = (Zi(0), . . . , Zm(0)) of particles 
withP (Z(0) ^ 0) > and having a random initial size $(0) of the final product. 
If 

ntVl . :^t, 



E 



for any t e (0, k), then 
if and only if x € (0, k). 



iiz(o)ir^^+$*(o) 

E$^ < 00 



< oo 
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Proof. Let be the total size of the final product produeed by all de- 
scendants of the k-th particle of type i of the zero generation. Clearly, the 
accumulated amount $ of the final product can be written as 

i=l k=l 

Given Z(0) = z = (zi, . . . , Zm) we have for a; < 1 

m Zi m 

E[$^|Z(0)=z] < E$^(0) + ^^E$?fc =E$^(0) + ^2iE[$^|Z(0) =ei] 

i=i k=i 1=1 
< E$^(0) + ||z|| max E[$^|Z(0) =ei], 

l<i<m 

while for a; > 1 there exists a constant R* such that (see, for instance. Theo- 
rem 5.2, page 22 in [27]) 



E[$"|Z(0)=z] < 

< (m+lf 

< (m + 1)^ 



E$^(0) + ^E^$,J 
i=i \fe=i / 

m 

E$^(0) + [^""IZW = Bi] 



i=l 



E$-(0) + R* ||z|r m.^ E [$-|Z(0) = e,] 

l<2<m 



By the total probability formula and the estimates above we obtain for x G 
(0,k): 

E^'^ = P(Z(0) =0)E$^(0)+ P(Z(0) =z)E[$="|Z(0) =z] 

zez+\{0} 



< ii'2E$^(0) + ii'3 max E[$^|Z(0) =ei]E ||Z(0)|| 

l<2<m L 



xVl 



< 00. 



On the other hand, for any z G N[)"\ {0} such that P (Z(0) = z) > 

> E [$"'|Z(0) = z] P (Z(0) = z) 
and the desired result for a; > k follows from Theorem 6. 



8 Polling systems with zero switchover times 
and MBPRCRE 

We start this section by the description of a connection between the BTPSFPRE 
with zero switchover times and the MBPRCRE. 

Consider a polling system with zero switchover times and assume that the 
operation of the initially idle system starts at the moment when a customer 
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arrives to a station J £ {1,2,..., m}. We would like to study the distribution 
of the busy period of the server which performs J — 1 switches of zero length 
and then starts the service of the customer arrived to station J. 

The length of this busy period is constituted by the time intervals spend by 
the server to make a random number of complete cycles (1^2— >1) 
plus the time needed to perform the last (may be, incomplete) route just before 
the moment when there are no customers in the system for the first time. 

Assume that in the course of the n— th service cycle within the initial busy 
period the service discipline at station i satisfies the branching property with 
vector-valued m.p.g.f. 



.^„(s;A) = (</.«(s;A),...,<^(r)(s;A)) 



where 



)(s;A) :=E 



1,2,. ..,m. 



Suppose that the sequence 0o(s;A), (/>i(s;A), ... is selected at random in an iid 
manner. Let further, for each customer, say j, served at station i during the 
n— th route the final product 4>i{n,j) and the vector {9ii{n,j), . . . , 0im{n,j)) of 
the numbers of new customers arrived to the system during the service time 
Ti{n,j) of the customer under consideration have the property 

{Oii{n,j), Oim{n,j); (j)i{n,j)) = (^ii(n), . . . , dim{n); (l)i{n)) . 

Here the final product may be not only Ti{n,j) but any nonnegative random vari- 
able being either dependent on the the tuple {9ii{n, j), . . . , Oim(n,j); Ti{n,j)) or 
independent on the performance of the system at all. Set h'i^} (s) := (j)n^ (s; 0), i = 
1, 2, . . . , m and for n = 0, 1, 2, . . . introduce m.p.g.f. 's 



F«(s;A)=E 



. g?im(n)g-Avi(n) 



and p.g.f.'s 



„Cii(")„«.2(n) 



' -s. 



by the equalities Fji^\s; A) = (^l'"^(s;A), 

F«(s;A) = # (s; A) , . , 

and/^)(s) = /i(r)(s). 



,Fi"')(s;A);A),z<m, (50) 



•,/^^ (s)) 



i <m. 



(51) 



We would like to describe conditions on the branching type polling system 
under which power moments of the amount of the final product accumulated in 
the system during its busy period are finite or infinite. Our results are based on 
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the following a bit long but important statement revealing connections between 
the behavior of certain characteristics of the busy periods of BTPSFPRE and 
related characteristics of MBPRCRE whose population and the size of the final 
product at moment are random. 

Theorem 17 The joint distribution of the number of customers at different 
stations at the end of the n—th service cycle and the amount of the final product 
accumulated in the system to the end of the n—th service cycle in a BTPSFPRE 
starting by a single customer at station J and stopped at the end of the first 
busy period coincides with the joint distribution of the number of particles in 
the n—th generation and the total amount of the final product produced for the 
n — 1 generations in a MBPRCRE whose generation is specified by a random 
number of particles and a random amount of the final product with m.p.g.f. 
Fq'^''(s;A) and where the joint distribution of the number of direct descendants 
and the amount of the final product produced by particles of different types of 
the k-th generation is given by the vector-valued m.p.g.f. 's 

Ffe(s;A) = (f1:'\s;X), Fi™^(s;A)) , k = 1,2, . . . ,n. 

To prove this theorem one should repeat almost literally the proof of Theo- 
rem 4 in [36] and we omit the respective arguments. 

Note that if, instead of a single individual at station J we would initially have 
in the system a batch of customers (fci, . . . , k„i) with ki customers at station i 
then the distribution of the initial number of particles and the initial size of 
the final product in the corresponding MBPRCRE should be specified by the 
m.p.g.f. 

m ^ 

Fo(s;A)=n(^o''M)) 
j=i 

We call the MBPRCRE described by Theorem 17 the associated MBPRCRE 
for the BTPSFPRE. 

Now we may reformulate the results of Section 3 in terms of our polling 

system. 

Let An := be the matrix with elements 

aij (n) := ^/^'^ (s) |s=i = Ef (n) (52) 
and let if„ := {hij{n))^.^^ be the matrix 

hij{n) := ^/li'Hs) |s=i = Eh%(n). 
Then in view of (51) amj{n) = hmj{n),j = 1, 2, . . . , m and for i < m 

aij{n) = hij{n)I {j < i} + ^ hik{n)akj{n). (53) 

k=i+l 
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For i = 1, . . . , m introduce auxiliary matrices 



1 




hii{n) 










i)(n) hiiin) /ii(i+i)(n) 













him iP') 







1 



where for each i the elements of the matrix _ff„ are located only in row i of the 
matrix Hn\ It is not difficult to check by (53) that 

= (54) 
Further, let C„ := (Ci (n), . . . , Cm(n))' be a random vector with components 



Ci{n) := j^F^\-^\^) |a=o = ^y^M), i 



1, . . . ,m 



and let c„ := (ci(n), . . . , c„i(n))' be a random vector with components 

Ci{n) := ^(/)W(l;A) \x=o = ^4>Mn), « = l,...,m. 
Then, by (50) Cm{n) = Cm{n) and, for i < m 



Cj(n) = Cj(n) 



fc=s+l 



hik{n)Ck{n). 



Hence we get C„ = c„ + -ff^C„ where 

F;^:=(/i„(n)/(z<j))r,=i 

is the upper triangular matrix generated by Thus, C„ = (iJ — iJ^) ^ c„. 

The next two statements are easy consequences of Theorem 16 and 17. 

The first theorem gives conditions under which a busy period of a BTPSF- 
PRE is infinite with positive probability. 

Theorem 18 Assume that the MBPRCRE associated with a BTPSFPRE is 
such that its underlying MBPRE satisfies conditions of Theorem 5 with a > 
and, in addition, condition (16) is valid. If<^ is the total size of the final product 
accumulated in the BTPSFPRE during a busy period then P (4> = oo) > 0. In 
particular, if the service time of any customer at any station is positive with 
probability 1 then the busy period of the BTPSFPRE is infinite with positive 
probability. 



30 



The statement of the results for a BTPSFPRE whose associated MBPRCRE 
is subcritical requires more efforts. 

Theorem 19 Assume that the MBPRCRE associated with a BTPSFPRE is 

subcritical, satisfies conditions of Theorem 6 and 

min Fy\O;0)>0. 

1< J<m 

If the parameter k specified by (20) is such that 



max E 

l<J<m 



< oo 



for any t G (0,«;), then there exists a constant C G (0,oo) such that 

P > y) ~ Cy-~, y ^ 00. (55) 

In particular, if the final product of any customer is its service time then the 
tail distribution of the length ^ of a busy period of the system satisfies (55). 

Corollary 20 Under the conditions of Theorem 19 E$^ < oo if and only if 
X G {0,k). 

Let us come back to Examples 1 and 2 considered at the beginning of the 
paper. 

Differentiating (2) at point s = 1 we see that the matrix if„ and the vector 
C„ in Example 1 have elements 

hij{n) = jij{n) + e^j{n)E [Ti(n)|Tj„] ,i,j = l,...,m 

and 

C,(n) = E[ri(n)|Ti„] ,i = l,...,m, 



while by differentiating (5) at point s = 1 and taking into account (3) we 
conclude after evident transformations that hu{n) = 0, i = l,...,m and, for 

, _ lain) (1 - 7m W) + Eij(n)E [Ti(n)|T^„] 

l-7u(n)-£,,(n)E[r,(n)|T,„] ' ^'^^ 

and hij{n) = oo, otherwise. 

Note that if the quantities eij{n), i,j = l,...,m are nonrandom, 7ii(n) = 
with probability 1 for alH = 1, . . . , m, and Tj„ is an exponential distribution 
with random parameter relations (56) and (57) look as follows: 



Hjyt'J — — ^ 



hij{n) 
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and 

^J■in - eii{n) > 0. 

They are in complete agreement with the respective formulas and restrictions 
of Sections 1.1 and 1.2 in [34]. 

Concluding remarks. Our results give a criterion allowing to answer the 
question: when the total amount of the final product accumulated in the 
polling system with zero switchover times during a busy period has finite or 
infinite moment of order x7 In this respect Theorem 19 refines and extends in 
several directions Theorem 1.1 in [34]. For instance, we do not require exponen- 
tiality of the service time distributions of customers and prove the mentioned 
criterion for a wide class of polling systems which arc not covered by the results 
of [34]. Moreover, we even describe the behavior of the tail distribution of 
Unfortunately, such a refinement is achieved for the expense of transparency of 
the conditions involved. The most essential of our hypotheses is Condition T, 
whose validity is established up to now only for a restricted class of nonnegative 
random matrices. The extension of the class of nonnegative random matrices 
for which Condition T is valid is an interesting and challenging problem. Some 
statements related with such circle of problems have been obtained quite re- 
cently in a number of papers (see, for instance, [17] and [26]). Unfortunately, 
they do not fit the case of measures concentrated on nonnegative matrices only. 
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